# t 


ICASE 



A SPLINE-BASED PARAMETER AND STATE ESTIMATION TECHNIQUE 
FOR STATIC MODELS OF ELASTIC SURFACES 


//O ^ 4 ^ 


H. T. Banks 
P. L. Daniel 
E. S. Arastrong 


Report No. 83-25 
June 28, 1983 


COPi 

NOV 2 0 

UANGLC. research CENTER 
UBRARV. NASA 

HAMPTON- 


r 

T — 

(NASA-CR- 180449) A SPLINE-BASED PAEANBTEE 

H87- 30107 



AND STATE BSTIHATION TECHNIQUE FOR 

STATIC 




MODELS OF ELASTIC SURFACES (NASA) 

62 P 




Avail: NTIS HC A04/MF A01 

CSCL I2A 

Unci as 


V 


G3/64 

0064388 

* 


NASA Langley Research Center, Hampton, Virginia 23665 


Operated by the 


UNIVERSITIES SPACE 



RESEARCH ASSOCIATION 



A Spline-Based Parameter and State Estimation Technique 
For Static Models of Elastic Surfaces 


May, 1983 


H. T. Banks 

Lefschetz Center for Dynamical Systems 
Division of Applied Mathematics 
Brown University 
Providence, RI 02912 



P. L. Daniel 

Department of Mathematics 
Southern Methodist University 
Dallas, TX 75275 



E. S. Armstrong 
Spacecraft Control Branch 
Flight Dynamics and Control Division 
NASA Langley Research Center 
Hampton, VA 23665 


Abstract 

We discuss parameter and state estimation techniques for an elliptic 
system arising in a developmental model for the antenna surface in the Maypole 
Hoop/Column antenna. A computational algorithm based on spline approximations 
for the state and elastic parameters is given and numerical results obtained 
using this algorithm are summarized. 
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I. Introduction 


Proposed large space structures have many characteristics that make them 
difficult to analyze and control [1]. They are typically highly flexible-- 
with components mathematically modeled by partial differential equations or 
very large systems of ordinary differential equations. They have many resonant 
frequencies-'possible low and closely spaced. Natural damping may be low and/or 
improperly modeled [2]. Coupled with stringent operational requirements of 
orientation, shape control, and vibration suppression, and the inability to 
perform adequate ground testing, these characteristics present an unconventional 
control design problem to the systems theorist. 

Frequently, linear multivariable control theory is applied to reduced 
order models obtained from finite element software with the goal of controlling 
the structure so as to meet the operational requirements without overly exciting 
higher order modes. Difficulties can arise essentially from the process of 
analyzing and designing controllers for flexible, inherently distributed systems 
employing lumped parameter models and finite dimensional control methodology. 

The effects of spillover, modeling error, and insufficient structural damping 
are well known [1, 3] and may force the designer into low authority control 
laws with compromised operational performance. 

An alternative approach would be to acknowledge the distributed nature 
of the problem and design the control system using distributed parameter 
techniques. Many aspects of distributed parameter control theory analogous 
to multivariable theory exist or are currently being developed [1, 4-7, 7a, 7b, 
7c]. Unfortunately, control laws derived from distributed parameter theory 
usually are infinite-dimensional and often require significant simplification 
to be realized with current sensor and actuator technology. Although distributed 
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parameter control methodology may ultimately be applied to large space structures, 
it is doubtful that first generation I'^'-ge space structures will have sufficient 
hardware and computer capability to apply the resulting solutions. 

However, even at the present time, much can be gained from the distributed 
parameter approach in the area of parameter estimation. Whether the lumped 
or distributed parameter control design approach is taken, there remains the 
underlying need for design models which adequately reflect the structural 
characteristics and which can be quickly and easily modified to carry out 
parametric studies [3] and to perform sensitivity and robustness analyses 
essential to control design [8]. There are also open questions in the under- 
standing and modeling of damping [2, 3] which may be better treated within a 
distributed parameter formulation. Finally, a distributed model appears to have 
the potential of facilitating the parameter estimation problem for large scale 
systems since quantities to be estimated will usually appear explicitly within 
coefficient functions of the partial differential equations. 

This report summarizes some of the results from an ongoing Langley Research 
Center program directed towards developing parameter estimation techniques for 
flexible systems modeled by partial differential equations with an emphasis on 
large space structures. The intent of the program is to produce general purpose 
techniques with a sound theoretical basis which are computationally efficient 
while contributing to Langley's technology development program in large space 
antennas [9]. Of the many techniques available for parameter estimation in 
distributed systems (for example, see [lOj) trie spline based estimation techniques 
of [11-15, 15a] appear well suited for large space structures applications and 
are currently being developed to treat this class of problems. Simultaneously, 
an estimation problem associated with the Maypole (Hoop/Column) antenna [16] 
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is being formulated and will be solved as part of the developmental process. 

The next section of this paper describes the Hoop/Column antenna and presents 
the identification problem being considered. The parameter estimation approach 
is then outlined and discussed in the context of the Hbop/Column application. 
Subsequent sections include mathematical details of the antenna application 
and numerical results. 

II. The Maypole (Hoop/Column) Antenna 

One of the planned activities of the NASA's Space Transportation System 
is the placement in earth orbit of a variety of large space antennas. Potential 
large space missions for the next two decades will require antennas and structures 
ranging from 30m to 20km in size. Applications include communications (mobile, 
trunking, etc.), remote sensing (soil moisture, salinity, etc.), deep space 
network (orbital relays), astronomy (x-ray observatory, optical array, radio 
telescope, very long baseline interferometry, etc.), energy, and space platforms. 

For the purpose of technology development, the NASA Large Scale Systems 
Technology (LSST) program office has pinpointed focus missions and identified 
future requirements for large space antennas for communications, earth 
sensing, and radio astronomy [9]. In this study, particular emphasis is placed 
on mesh deployable antennas in the 50-120 meter diameter category. Communication 
satellites of this size will require a pointing accuracy of 0.035° and surface 
accuracy of 4-8min. One such antenna is the Maypole (Hoop/Column) antenna shown 
for the 100m point-design in Figures 1, 2, and 6. This antenna concept . 

has been selected by the LSST office for development by the Harris Corporation, 
Melbourne, Florida, under contract to the Langley Research Center [16]. 
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The Hoop/Column antenna consists of a knitted gold-plated molybdenum 
wire reflective mesh stretched over a collapsible hoop that supplies the 
rigidity necessary to maintain a circular outer shape. The mesh grid can be 
varied to meet a given radio frequency reflectivity requirement. The annular 
membrane-like reflector surface surrounds a telescoping mast which provides 
anchoring locations for the mesh center section (Figure 2). The mast also 
provides anchoring for cables that connect the top end of the mast to the outer 
hoop and the bottom end of the mast to 48 equally spaced radial graphite cord 
truss systems woven through the mesh surface [16]. Tensions on the upper 
(quartz) cables and outer lower (graphite epoxy) cables are counter balanced 
to provide stiffness to the hoop structure. The inner lower cables produce, 
through the truss systems, distributed surface loading to control the shape 
of four circular reflective dishes (Figures 2 and 6) on the mesh surface. 

Flat, conical, parabolic, or spherical dish surfaces can be produced using this 
cable drawing technique. 

After deployment or after a long period of operation, the reflector surface 
may require adjustment. Optical sensors are to be located on the upper mast 
which measure angles of retroreflective targets placed on the truss radial 
cord edges on the antenna surface. This information can then be processed 
using a ground-based computer to determine a data set of values of mesh surface 
location at selected target points. If necessary, a new set of shaping (control) 
cord tensions can be fed back to the antenna for adjustment. 

It is desirable to have an identification procedure which allows one to 
estimate the antenna mesh shape at arbitrary surface points and the distributed 
loading from data set observations. It can also be anticipated that environ- 
mental stresses and the effects of aging will alter the mesh material properties. 
The identification procedure must also allow one to address this issue. 
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It is the intention of the authors to develop identification procedures 
for use in distributed parameter models of the mesh surface. Considering the 
antenna to be fully deployed and in static equilibrium, we are currently 
deriving a mathematical model which describes the antenna surface deviation 
from a curved equilibrium configuration (for preliminary findings, see [16a]). 
Using a cylindrical coordinate system with the z axis along the mast, it is 
expected that the resulting model will entail a system of coupled second order 
linear partial differential equations in two spatial variables. The coefficients 
of these equations are functions of the material properties of the stretched 
mesh. The derivation and computer software for this model are still under 
development. In the meantime, a simpler developmental (prototype) problem 
has been solved which is descriptive of the original problem and for which the 
software produced will hopefully be readily extended for use in the more general 
case. 

For the developmental problem, the loading is assumed to be normal to the 
plane containing the hoop rim and the mesh surface is assumed to be described 
by the static two-dimensional stretched membrane equation [17] with variable 
stiffness (elastic) coefficients and appropriate boundary conditions for the 
Hoop/Column geometry. (We note that there is ample precedent in the literature 
[17a, 17b] for use of a scalar Poisson's equation for the two-dimensional 
stretched membrane as a prototype or developmental model in the study of large 
space antennas.) Mathematically, in polar coordinates, we have 


1 A. 

r 3r 




3U- 


( 1 ) 


where u(r,e) is the vertical displacement of the mesh from the hoop plane. 


6 


f(r,e) is the distributed loading force per unit area, and E(r,e) > 0 is the 
distributed stiffness (elastic) coefficient of the mesh surface {force/unit 
length). Equation (1) is to be solved over the annular region n = [e,R] x- 
Appropriate boundary conditions are 

u(c,e) = Uq 

u(R,e) = 0 ( 2 ) 

u(r,0) = u(r,2ir), 

where R is the radius from the mast center to the circular outer hoop, e is the 
radius from the mast to the beginning of the mesh surface (see Figure 6), and Uq 
is the coordinate at r = e of the mesh surface below the outer hoop plane. For the 
100-meter point design antenna, one finds 

R = 50 m 

e = 8.235rn 

Uq = - 7. 5 m . 

We further assume that the distributed loading along with a data set of 

vertical displacements, at selected points (r.,e.) on the mesh. 

* 

surface is known. Given this information, the developmental problem is to 
estimate the material properties of the mesh as represented by E(r,e) and 
produce estimates of the surface represented by u(r,e) at arbitrary (r,e) 
points within n. The procedure applied to solve this problem is discussed 
in the next section. 
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III. The Parameter • Estimation Approximation Scheme 

The first two authors and thei'" colleagues have derived techniques for 
approximating the solutions to systems identification and control problems 
involving delay equation models and partial differential equation models in 
one spatial variable and have used them in a variety of applications [18, 18a]. 
The Hoop/Column application requires an extension of the theory and numerical 
algorithms to el 1 iptic distributed systems in several spatial variables. The 
approach, when specialized to the parameter estimation problem, may be summarized 
as follows. (1) Select a distributed parameter formulation containing unknown 
parameters for a specific system. (2) Mathematically "project" the formulation 
down onto a finite dimensional subspace through some approximation procedure 
such as finite differences, finite elements, etc. (3) Solve the parameter 
estimation problem within the finite dimensional subspace obtaining a parameter 
estimate dependent upon the order of the approximation embodied in the subspace. 
(4) Successively increase the order of the approximation and, in each case, 
solve the parameter estimation problem so as to construct a sequence of parameter 
estimates ordered with increasing refinement of the approximation scheme. (5) 
Seek a mathematical theory which provides conditions under which the sequence 
of approximate solutions approaches the distributed solution as the subspace 
dimension increases with a convergent underlying sequence of parameter estimates. 

In applying this approach to the developmental problem described in this 
report, the stiffness function is parametrized in terms of cubic splines 
which converts the estimation of E(r,e) into a finite dimensional parameter 
estimation problem. After writing the energy functional generic to the 
membrane equation, we use the Galerkin procedure to project the distributed 
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formulation onto a finite dimensional state subspace spanned by tensor products 
of linear spline functions defined over Q. The approximate displacement thus 
obtained is expressible in terms of the spline basis functions. The Galerkin 
procedure in this case yields algebraic equations which define the displacement 
approximation coordinates in terms of the unknown parameters. In order to 
solve the approximating parameter estimation problem, the parameters defining 
E(r,e) are chosen so that a least squares measure of the fit error between the 
observed and predicted (by the approximate state) data set is minimized. Finally, 
following steps (4) and (5) an algorithm is constructed to determine the 
order of the linear spline approximation above which little or no further 
improvement is obtained in the unknown parameters as one increases the dimension 
of the subspaces. 

Following this procedure, a one-dimensional version of the developmental 
problem has been solved [19] in which the schemes proposed were successfully 
tested. Further details of the two-dimensional case follow. 

IV. Finite Dimensional Approximations 

We choose a Galerkin procedure [20, 20a] with linear spline basis functions 
to perform the finite dimensional approximation for the developmental problem 
in which E(r,e) and u(r,e) of equation (1) are to be estimated. The boundary 
conditions (2) are first converted to homogeneous form by introducing the new 
dependent variable 

y(r,e) = u(r,f.) - ( ) Ug (4) 


whereby equation (1) becomes 
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CrE(r,6) ^ [E(r.o) 5^] = f(r.o) ♦ i :^ [ 


r 3r 


9r 


90 - 


r ar 


rE(r,8)Uj 
e - R “ 


(5) 


with boundary conditions 
y(e,e) = 0 

y(R.e) = 0 (6) 

y(r,0) = y(r, 2ir). 

Following the standard formulation (see [20, 20a]) for the weak or 
variational form of (5), the energy functional E associated with (5) is 


E(2) 


= I 

0 


2n R 


!{i 


E(r,e)7z • vz - ?(r,0)z^ rdrdG, 


(7) 


where v is the gradient in polar coordinates which, in the form used here , 
is equivalent to 


/ 5 3 \T 

Dr ’ r 30 ' 


The function f is given by 


( 8 ) 


f{r,0) = f(r, 0 ) + ^ 


Dr 


rE(r,0)u 

c - R 


0 


(9) 


and the vertical displacement z(r,0) of the mesh surface away from the hoop 
equilibrium plane is a function satisfying the boundary conditions (6) and 
possessing first derivatives on o in the mean square sense (we denote this by 
z e Ro’per^*^^ ^ first variation 6E of E about the function y(r,0) 

is given by 


rdrdo 


2it R 


6£(y;v) = / / / E(r,o)vy • vv - ?(r,0)vl 


2ir . ^ C 'v 

= / / sE(r,8)vy • w - [f(r,e)v +£(r,o)k • 7v] ) 

0 c ^ J 


rdrde 



( 10 ) 


and V is an arbitrary function in Z = ov(^)- 

U» k“' 

Under sufficient smoothness assumptions one can integrate by parts in 
the first term in (10). Then, if one invokes the conditions that y provides a 
stationary value for E, the fact that 6E(y;v) vanishes for all v in Z yields 
equation (5). However, we shall not follow this course, but rather we shall 
apply the Galerkin procedure directly to the weak form of (5) given by 
6E(y;v) = 0. 

Given a finite dimensional subspace Z of Z, the Galerkin procedure defines 

the approximation y as the solution in Z of 

2irRr 2-n R f ^ 

/ / { E(r,o)vy • vv ) rdrdo = / f I f(r,o)v +E(r,e)k • vv 

0 e I j 0 t L 

for al 1 v € Z. 

For computational efficiency, the basis functions used for the representations 

of y in (12) are taken as tensor products of linear B-splines [20, p. 27; 20a, 

^ ' M N 

p. 100]. Thus V and y are in the space spanned by v^j , where 


^ rdrde (12) 
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v*?j'^(r,0) = a'J''{r)8j(o), (i = 1 , . . . , M-1 ; j = l,...,N), (13) 

where and have the following form: On the interval [c,R] (subdivided by 

defining partition points r*^ = e + (R-e)(j-l )/M, j = 1, .... Nh-l), a*!* is the 
i^^ linear B-spline basis element, i = 1 ,2, . . . ,M-1 , given in the following graphical 
representati on 



0 ^ (j-1) ^ j Ip ( j+1 ) 2it 


and define 6|^(e) as follows. 



0 ^ (;M)2t, 

N N 


12 


For y(r,e) within the subspace spanned by we can write 

J 

i=l j=l ’ ^ 


( 14 ) 


Note that the coordinate w^’.*^ is just the value of y^’^(r,e) at 

* J 


r = £ + 


i(R-e) 


0 = 

N • 


Replacing y(r, o) in ( 12 ) by y^’^(r,e) from ( 14 ), (and letting v(r,e) = v^j^(r,o) 
for i = 1 , M- 1 ; j = 1 , N) we generate the following set of algebraic 

M iN 

equations for w. j . 


j .^1 ''' ij ^(',ij 


(k= 1 , .... M- 1 ; r= 1 , N), ( 15 ) 


where 

2 Ti R 

C.ij "II E(r. 6 )v[a'J(r)f!j^e)]. 7 [t, 5 J(r)B^( 6 )]rdrde 

2tt r 

= ^ / E(r,e) 6 j( 8 )Bj( 6 )[ <.”(r)][ |:a”(i-)]rdrdo ( 16 ) 

+ / / E(i-.o) -J— j.J' — [ 6 j( 8 )][ ^ B”(e)]rdrde . 

M N 2 Ti R » J 

= / / f(r,0)(i|^(r) Bp(G)rdrdo, (17) 


and 
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2tt R 

~ / E(r,e)k6'^(e)[ ^ a^(r)]rdrde. 


We note here that 


.MN MN MN 

•^k^.ij ■ Ne,kj ■ ■ Nj.kf • 


The systen (15) can be written in vector matrix notation as 


.MN .MN ^MN 

A W = C 


where, in partitioned matrix form. 


= (Ap^), (p -- 1, ..., M-1; q = 1, ..., M-1) 


and, for fixed p and q. 


a"N = (k”N .) 
pq P'.qj 


for (row) i = 1, ..., N, and (column) j = 1, ..., N. From properties (19), 
we also have that 


. [A*f . a«N 

pq pq-' qp 


MN MN 

In addition, the vectors w and c are given by 
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For large values of the product MN. is high order, ( (M-1 )N x (M-1 )N) . 
It is a symmetric, sparse matrix which can \v partitioned as (M-1)^ NxN 
submatrices. Sparse matrix techniques can be avoided in solving (20) if we 
assume that E(r,o) is separable. That is, 
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E(r,0) = E^(r) £ 2 ( 0 ) (26) 


where E^(r) 

> 0 and £ 2 ( 0 ) > 0; we may now write 



„MN 

V, 

. = a", + c^. d"" 

.qj ji qp ji qp 

(27) 

where 

all. 

01 

2tt 

= a^.. = / E2(0)0'J(0)Bl'(0)d0 , 

(28) 


qp 

' C “ ( ^ ^ »p(r)]rdr . 

(29) 


N 

^ji 

= - l’ E2(»)[ i e>)][ ® 8?(e)3<‘» 

(30) 

and 

d"^ 

qp 

^ R ot^(r) Op(r) 

£,(-■) dr. 

e 

(31) 


Utilizing (27) and defining matrices 


'V'M 

= 


= 

1» •••* Nj j ~ Ij •••> 

N),. 

(32) 

'^M 
B = 

<c> 

(q = 

1 5 • • • > M” Ij p-*lj •• 

M-T) , 

(33) 

'^'N 

= 

<^u' 

(i = 

1, j"l> ••• 

, N) . 

(34) 

= 

(d” ) 
' pq' 

(P = 

1, .... M-1; q = 1. .. 

.. M-1), 

(35) 

we are able 

to rewrite 

(21) 

and (22) as 
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pq qp pq 


p !^ 


where the symbol © denotes matrix direct product [21 ]• 
Equation (20) now becomes 

which is equivalent [21, p. 261] to the matrix equation 
wMN ^ ^MN ^ |MN 


(36) 

(37) 


(38) 


(39) 


where 


= (w"J) 


(i = 1, .... M-1; j = 1, N) 


(40) 


and 


^MN 


= (pMN 

' 1J 


+ G^^) 

ij 


(i = 1, 


M-1; j = 1. 


N) 


(41) 


MN 

The coefficient matrices of W in (39) have numerically attractive properties 
(i) all are symmetric, (ii) B , C , and D are banded (tri-diagonal), and 

, . . . , 'V.M -VM 

(lii) A and D are positive definite. 

Research to construct a numerical algorithm for solving (39) which 
utilizes these properties is planned. At present, (39) is rewritten in the 
equivalent form 






] = (d”) 




( 42 ) 


and solved by the Bartels-Stewart algorithm [22]. 


18 


In order to estimate, via a numerical schene, the functional coefficients 
and £ 2 * we must further parametrize these functions so that identification 
is performed over a finite-dimensional (instead of an infinite-dimensional) 
parameter set. Our approach here is to consider a spline-based representation 
for each so that the shape of E^, E 2 is not assumed in advance (as would be 
the case if E-|, E 2 were assumed to be affine functions, for example). To this 
end, we let 

= I (^ 3 ) 

' k=l 

Ep(o) ^ I i5.m-((!) (44) 

^ j = l J J 

where and 6j are scalar parameters and and are cubic B-spline functions 
defined [20, p. 61] over [e,R] and [0,2tt], respectively, whose orders are 
independent of M and N. The functions have been modified so that y. satisfies 
periodic boundary conditions. Using (43) and (44), 


= I 6. A' 


k=l 

M 


k "k 


1 'V', 


B" = I V, B' 


k=l 

^1 


M 
k "k 


r = 




k -1 


M 


1 'V'. 


D = I v.D' 


k=l 


k\ 


(45) 


(46) 


(47) 


(48) 
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where the following matrices have been defined: 


^ U|^(e) 6 p( 0 ) 6 q{ 0 )de) ) (p=l, .... N; q = 1, .... N) (49) 




(50) 


(p=l, M-1; q = 1, .... M-1) 


%i / 2n 




(51) 


(p = 1 , . . . , N; q = 1 , . . . , N) 


M, . M, V 


-(l A,(r) 


(52) 


(P= 1. 


M-1; q = 1 M-1). 


We also obtain, using (43) and (44) to characterize representations for and 



(53) 


M-1; j = 1 


N). 
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We turn next to implementation of the parameter estimation scheme, i.e., 
the numerical determination of k = 1, . . . , , and 6^., j = 1 , N^, 

that appear in (43) and (44) and correspond to "optimal" values of the parameters 
E| and £ 2 * 

V. The Parameter Estimation Algorithm 

Appealing to the ideas found in previous sections, we may now detail an 
algorithm for estimating the coefficients V|^, k = 1, ..., M^and 
3 = 1, ..., N^,for E(r,o) that provide the "best fit" between approximations 
for the state u and observed data u^ obtained from various sample points on 
the surface. We may equivalently consider data for y by making the trans- 
formation 

for i = 1, . . . , L^, and j = 1 , . . . , Lq. 

We organize the parameter estimation algorithm into the following steps. 

1. Select an order of approximation for the cubic spline elements Aj^, 

k = 1, ..., M.J , and j = 1, ..., N.j, used to represent E^ and E 2 . 

Set n = 1 . 

2. Select M and N, the orders of the linear spline basis elements used 
to represent u^’^ (and 

3. Assume a nominal set of values for 

V = (v-| , V2> . . . , V|^ ) . (55) 

and 


(56) 
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4. 


5. 


Calculate the coefficient matrices in (42) and solve for . 

Calculate, from (14), y^’*^(r .,e .;v,6) and evaluate 

* J 


^r ^0 
(v,6) =11 

i=l j=l 




( 57 ) 


6. Proceed to step 8 if J^’*^(v,6) is sufficiently small. Otherwise, 

through an optimization procedure, determine a new pair (v,^) which 
decreases the value of If no such pair can be found, go to 

step 8. 

f\j 

7. Set (v,5) = (v,6) and return to step 4. 

8. Preserve the current values of and the corresponding (v,6) pair 
as the n entry in a sequence of these pairs, ordered with increasing 
M and N. 

9. Proceed to step 10 if sufficient data has been obtained to analyze 
the sequences. Otherwise, set n = n+1 and return to step 2 with 
increased M and N. The current values of (v,6) will be used as 
initial values for the next optimization process. 

10. From analysis of the numerical sequences, select the M, N entry 

which indicates the best numerical results. The corresponding (v,6) 
pair yields E(r,e) which determines the material properties of the 
antenna mesh. The matrix W™^(v,5), when used in conjunction with 
(14), determines an approximation y^’^ of the shape of the antenna 
surface. 

For the examples reported in the next section, a Levenberg-Marquardt [23] 
(Damped Least Squares) nonlinear programming scheme was applied at step 6. 
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A convergence theory has been developed for our parameter estimation 
algorithm and will be outlined in section VII below. We instead focus next 
on a comprehensive discussion of our numerical findings to date. 


VI. Numerical Results 

We present here our findings for a number of test examples. In Example 1 

arid 2 we demonstrate the effectiveness of the scheme outlined in Section V 

as it is applied to representative model equations. We then turn to Example 3, 

which is constructed based on proposed physical characteristics of the antenna 

surface (such specifications are detailed in [16]). Lacking at present 

actual experimental data from the Hoop/Column antenna, we hope, in this 

example, to best emulate the type of data expected from the antenna. 

For each example, a L^x grid {y^^jCr. ,6^ ), i = 1, ..., L^, j = 1, ..., Lg} 

of the data points is generated by choosing a desired shape u and evaluating 

u (with noise added in some cases) at (r.,e), i = 1, ..., L , j = 1, ..., L ; 

1 J r 0 

the numerical package is then given the data and the distributed load f, where f 
is calculated by substituting the noise-free u and a choice of elastic parameter 
E into equation (1). For desired values of N and M, and an initial guess E^ for E, 
we apply standard optimization schemes to minimize (given in (57)). In 
each case, we take E(r,o) = E-|(r)E2(0) where we approximate E^ and E 2 by 
linear combinations of cubic spline functions in r and 0, respectively. 

Throughout, equations (43) and (44) are used to represent E-j and E 2 , using 


= 4 . 
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All numerical examples presented here were computed on the CDC 6600 at S 

Southern Methodist University. We would like to express our sincere appreciation 

to David Krakosky at S.M.U. who assi^ited in the preparation of data and 

graphical display of these results. The optimization scheme employed to 
• - • MN . 

minimize J is the IMSL version (ZXSSQ) of the Levenberg-Marquardt algorithm, 
where we typically use default values of IMSL parameters. 

For each example reported below, values of M (or N) are specified; in 
order that M- 1 basis elements are used in both r and 6 directions, we take 
N = M - 1, Two measures of performance will be given in each case: 

(which provides a good measure of state approximation) and where, using 
I*! to denote the L2 norm on [c,R]x [0,2ti], 


.M,N 


IE^N. E, 


100 % 


measures the relative error between the "true" parameter E and the "optimal" 

MN 

parameter E associated with the (M,N)th approximate parameter estimation 
problem. We also report CP time which can be reduced considerably by relaxing 
the convergence criteria for the IMSL package and by taking as initial guesses 

'VA/ 

(for the (M,N)th problem) the resultant E^^^ from an earlier run (M < M, fj < N). 
This procedure is used in Example 3. 

Example 1.1 : Here we take e = 1, R = 6 , Uq=- 50, u(r,e) = 

-(6-r)[(r-l)(sin(2e+2) + 10] and E(r,e) = .25(r-3)^ +2, for 1 < r < 6 
and 0 £ 0 £ 2ti. A 7 x 7 grid of sample data is computed; throughout E2 is 
fixed at E2 H 1 so that we identify E-j only, starting from an initial guess 
of E° = .001. Our findings are summarized in Table 1 and Figure 3. 
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Example 1.2 : We repeat the same example but now add random noise (distributed 

normally with mean 0 and variance .75) to the data, resulting in approximately 
5% relative error. Numerical and graphical results are found in Table 2 and 
Figure 4* 

Example 1.3 : We repeat Example 1.1 (i.e., noise-free data) except that now 

we estimate E 2 ( 0 ) as well. We initially set E° = .1 and E 2 (e) = . 359(2 - sin(-|) ) 
and hold fixed the leading coefficient (in the spline expansion (44)) for 
E 2 (o). The results for M = 25 are found in Table 3 and Figures 5(a) and (b). 

Example 2 : In this example, we choose Uq=- 24, e = 2, R = 6 , 

u(r,e) = -(6-r)[(r-2)(cos2e+4)+6] and E(r,e) = (-^r + 7)(2 --^sin ^e), for 

r e [2,6] and e € [0,2n]. Again a 7x 7 grid of data is used in the optimization 

procedure. As an initial guess for E(r,o) = E^(r)E 2 (e). we let E°(r)E 2 ( 6 ) = 1 

and search for both E.j and E 2 , holding the first coefficient fixed in the 

cubic spline representation for E 2 * A summary of our findings for Example 2 

may be found in Table 4 below. 


EX AMPLE 1 


Table 1: Results for Example 1 


M 

jM,N 

(j*)' 


CP time (Sec.) 

4 

4.75% 

19.86 


8.52 

8 

19.64% 

3.61 


23.97 

16 

4.66% 

0.44 


89.43 

24 

2.25% 

0.23 


284.03 


32 


0 . 681 K 


0.11 


655.49 
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M 

Table 2: 

Results 

for Example 1 . 2 



CP time (Sec.) 

4 

21.83% 

21.26 

6.78 

8 

18.92% 

6.82 

20.80 

16 

4.44% 

5.73 

104.97 

24 

2.34% 

5.65 

410.51 

32 

1.17% 

5.64 

815.85 



Table 3: Results 

for Example 1.3 

M 

,M.N 

CP time (Sec.) 

25 

9.16% 0.17 

370.87 



Table 4: 

Results 

for Example 2 

M 

„M.N 

(jMN)% 

CP time (Sec.) 

4 

6.86% 

2.40 

11.03 

8 

3.41% 

1.58 

34.21 

16 

0.99% 

0.38 

313.68 

24 

0.48% 

0.10 

1200.00 
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Example 3 : We present here our efforts to construct functions u and E that 
follow many of the specifications of the Hoop/Column antenna surface found in 
[16]. 


As shown in Figure 6, the parent reflector has four separate areas of 
illumination or aperture on its surface. Each separate area is assumed to 
have the same parabolic shape given, for 0 £ e £ ^ and e £ r £ R, by 


Uf^(R-r) 

R - e 1 ] . 0 1 0 1 36 


u (r,9) = 


Uo(R-r) 


R - t; q](o) + 1] . 3'^6 - ® -- ^36 


(59) 


UnCR-r) , 

_0„ [k (-V) q3(0) . 1], 


17ir „ n 
36 - ® - 2 


where 


q^(e) = sine + cose 


(60) 


2 

q2(e) = ^2(6 - ^)^/6 + 2 (9 - (^) (®“ 36 ^ de^ * ^ 1 ^ 36 ^* 

de 


(61) 


- 279936 




( 0 ) 


.2 d'^qi 

2592 ,^2 


ddi ^ 


93 ( 0 ) 


®3^® 


17n)3/6,l (,.U.)2l^, 

d 


36 


17n. ^'’l ,17ir. 


'’l^ 36 ^ ’ 


(63) 



27 


279936 

_3 


{q,(f) - 


2592 




(Ujl\ jl ^ / 17ti x 
' 36 ^ ■ 36 de ^ 36 ' ■ 


(^)) 


n' 36 


(64) 


The parameter k > 0 is a stretch factor used to perturb the surface below the 
conic (k = 0) shape. 

For the complete surface, we define, for (- < r < R, 


ii(r,0) 


u°(r,o) 
u^(r,0 - j) , 
u^(r,.e - n) , 
u°(r,e - f-) , 


fi.<. 2 



0 < 71 


7T 




6 < 277 


(65) 


The cubic polynomial fits (61) and (63) are used to ensure smoothness in e, 
in regions near e = ^, n, 2n. 

We turn now to equations for E(r,e) = r^(r)E 2 ( 0 ). It is expected that 
the mesh will be stiffest near the outer hoop (r = R) and around the inner 
radius (r = c). For this reason we choose 

E^(r) = 2t - Tsin[n (c .1 .i R) (66) 

where t is a constant dependent on the mesh material. Stiffness in the angular 

/ 

direction is expected to be uniform, given here by 
£2(0) = 


T . 


(67) 


28 


From data provided in [16], it appears that a reasonable value for t (given 
in units /liTm ) is 

T = 3.391 ; (58) 

similarly, other parameters are estimated to be UQ=-7.5m, c = 8.235m, and 
R = 50m. For the example reported here, we take k = .25. 

For each example below, a 24x 24 grid of data points is calculated by 
determining u at (r^, 0^) where are equally spaced in [0,2tt], j = 1, .... 24, 
and the r^ are spaced in a non-uniform manner, with more data points concentrated 
near r = R. 

Example 3.1 : In an effort to determine how well the state approximation 

scheme performs when E is fixed at the "true" value E = E, we solved the 
"forward" problem initially to assess what values of to expect in later 
"inverse" problem runs. To this end, we compare in Table 5 below for 
various values of M (N = M-^). To interpret Table 5, it is helpful to calculate 

the average error at each of the 576 sampling locations: In the case of M = 4 we 

compute the average difference |u(r.,e.) - u*^’^(r. ,0 . ) | = {J^^*'*/576)'* = .056. 

Table 5: Example 3. 1 




4 

1.82 

8 

1.88 

16 

1.92 

24 

1.94 

30 

1.94 


We note that a little reflection leads to the conclusion that the values 
of J reported in Table 5 are not lower bounds for the values of J obtained in 



the minimization procedures detailed below (here we are using the approximate 
states and a spline approximation to the "true" E). This is evident in the 
subsequent tables for this example. 


Example 3.2 : We estimate E.j(r) caly, holding E2(e) fixed at the true value, 

£2(0)- For each M, the initial guess for E.| is E^(r) = 1. Our results are 
summarized in Table 6, and graphs comparing E^j*^ to E^ are found in Figure 7. 

Example 3.3 : We now hold E^(r) fixed, E^{r) 5 F.j(r), and estimate E2 from 

the starting guess of £2(0) = 1 + .5 cose (in the case of N = 4). For 
N = 8, 16, 24, and 30, we use the earlier converged value of E^^ as the 
initial guess (e.g., E^ for N = 8 is E^^^, N = 4). Findings for this example 
are summarized in Table 7; graphs of E^^, plotted against ^2, may be found 
in Figure 8. 

Example 3.4 : We estimate both E-|(r) and £2(0). Initial guesses for N = 4 

are given by E^(r) = 5 and £^(0) = 1 - | sine ; for N = 8, 16, 24, and 30, 
previous estimates (E*j*^, E^^) are used. In each case, the first coefficient 
for E2 is held fixed. The convergence results for this example may be found 
in Table 8. 

Example 3.5 : This example is identical to that of Example 3.2, except that 

A 

random noise has been added to the data. The new data u(r.,e.) is computed 

^ J 

by setting 

where n-- is a zero mean normally distributed random number with unit variance 

* \J 
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PAGE IS 

LXAMP]J 3 m POC« QUALmr 


N 

b«.n 


CP time (Sec.) 

4 

37.38‘-i 

0.904 

25.06 

8 

7.50^: 

0.679 

33. 73 

16 

4.06:1 

0.678 

168.08 

24 

4.64S 

0.668 

226.67 

30 

4.87^0 

0 . t> ^ b 

440.79 



L\.. 

3.3 

i\ 

,M,N 

• MNv’ 
) 

CP time (Sec.) 

4 

/. 74-; 

0.983 

14.37 

8 

6,20;. 

0.763 

19.31 

16 

6.244 

0. 763 

56.40 

24 

6.261. 

0.758 

73.68 

30 

6.254 

0.765 

133.07 



TaM e E v aii^ e 

3.4 

N 

,M,N 

(/.N,V 

CP time (Sec.) 

4 

43.274 

0.903 

26.88 

8 

5.674 

0.635 

127.31 

16 

5. 784 

0.636 

129.81 

24 

6.604 

0.627 

139.19 

30 

7.004 

0.634 

378.21 
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Thus, with probability of .95, the new data will satisfy 

|u(r.,e.) - u(r.,e.)| 

L_J_ < 05 

|u(r.,9j)i . 

so that data with, roughly, 5% noise level is used here. 

We take = 1.0 for N = 4 and 8, and previous estimates are used as 

initial guesses for N = 16, 24, and 30. Graphical representations of vs 
E-| are found in Figure 9; our findings are summarized in Table 9. 

Table 9: Example 3.5 

N (J^^) ^ CP time (Sec.) 

4 
8 
16 
24 
30 


148.86% 


5.05% 


6.76% 

4.99% 

5.33% 


2.669 

2.616 

2.612 

2.610 

2.611 


41.57 

46.73 

45.72 

169.22 


255.90 
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VII. Theoretical Convergence Results 

We present here a theoretical framework for the approximation scheme out- 
lined in Section IV. We shall first formulate the variational equation of state 
and investigate some of its properties; spline-based approximating equations 
(based on either linear or cubic spline elements) will then be constructed so 
that we may argue the convergence of approximate states and parameters to the 
"true" state and "optimal" value of the unknown parameter. 

Appealing to the notation of Section II, we may write the total energy of 
the antenna surface as 

E(u) = // (^ E |vul^ - uf)dsj 

'X, 

where the elastic coefficient E and distributed loading f are written in terms of 
Cartesian coordinates in n = {(x.j,X 2 ) 6R | < x^ + x^ < R } and 

u € = {u e H^(n) I u = y + U, ye Hq(o)}; here 

U(Xi,X2) = ((x^ + x^)’^ - R) , (x^.xg)^ fi . (69) 

From minimum energy considerations, v/e know that the state variable u is a 
stationary point of E. That is, E’(u;v) = 0 for all v € HQ(fi), where E* denotes 

the Frechet derivative of E. We thus write, for any ve Hq(?>), 

0 = E’(u;v) 

= // {Evu • vv - fvldfi , 

•X, 

n 

which yields the state equation for u Hg(?i), 
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afg(u,v) = <f,v>, V C Hj(ft) . ( 70 ) 

Here <*,*> is the usual L 2 (fi) inner product and 
a^(u,v) = // Evu • 7v d^. 

We remark that if u and E are sufficiently smooth, we actually may use Green's 
Theorem, ([27], p. 342) 

// EV4> ^ Vip = / E 4 > ds - // .tv (Evt) 

'h Mi ' "h 

(for t e h\ t e H^) to rewrite (70) as 

- // v7 • (Evu)do = // vfd?5 

since v€ Hq(o). Thus, for sufficiently smooth E and u, the state equation in u 
becomes Poisson's Equation with variable coefficient, 

- V • (Evu) = f on . ( 71 ) 

Before we consider various properties of the state equation (70), it is 
natural to first transform the state variable u into a new variable y satisfying 
homogeneous boundary conditions. To this end we define 

y = u - U 

and note that equation (70) may be rewritten in terms of y, v € Hq(o) as 


a^(y,v) = <f,v> - a^(U,v); 


(72) 
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further, if E € an integration by parts yields 


a^(y>v) = <F,v> (73) 

where F = f + v • (EvU). 

Throughout we shall assume that f€ L 2 (fi) and E € , where the set 

is a given set satisfying 

? {E€ I 0 < m < E(x) < m . |E|-, < W}. 

These conditions on E ensure that, for any y 6 Hglfi), 

'"|yli 5. a^(y.y) 

(1*1 . denotes the usual h 3(?2) norm) so that we may apply Theorem 6, p. 295 of [24], 
J 

to state the following existence theorem. 


Theorem 1 : There exists a unique solution u£ Hg(f$) to the generalized 

Dirichlet problem (70) (and thus there exists a unique solution y € Hj(fi) to 
(72)). 

In addition, we obtain the following regularity result, for E, f sufficiently 
smooth. 

Corollary 1: Let E € and f £ l 2 (^). The solution y to (72) satisfies 

y £ H^(^) and 

for some constant c«. If, in addition, f £- and E € C^, then y € and, 

for some constant c^ , 

|yl4 1 I2 + IEI3 + ly|) • 



35 


Proof ; Theorem 17.2, p. 67, of [25] may be applied to (73) to obtain 
lyl2+j 1 c(|p|j + |y|) 

1 c(|f|j + |v- (EvU)lj + iy|) 

for j = 0, 2, and some constant c, from which (75) and (76) follow. 

Clearly, statements analogous to (75), (76) may be made for the solution u 
to (70) since u y + U. In fact, although all results reported below involve 
y, it is easy to see that the same findings hold for u. We first establish some 
a priori bounds for y. 

Lemma 1 : Let E£ , f €. bounded in L2(f0* 3nd let y be the solution 

to (72). There exist constants Kq, (depending on and only) such that 

|y(E.n| 1 Kq (77) 

and 

|vy(E,f)l<K^ (78) 

uniformly in E e Q.| , f €, 

Proof : For v £ hJ and y = y(E,f), 

l'a[:(y.v)| . |F||v| 
so that, letting y = v, 
m |vy| < a^(y,y) 


1 |F||y| 
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where [Fj is bounded uniformly in E € and f € F-j. It thus follows that 

for some constant kg and, applying Poincare's inequality (given by 
ly|<k^lvy| 

for some constant k-j and y € Hq*, see for example, pp. 158-159 of [26]) we obtain 
(78), 

|vy| < k^k^ = K^. 

Combining (80) and (78), inequality (77) also obtains, where Kq and are 
determined independently of f€ F^, E£ Q^. 

We turn now to the problem of approximating the state variable y for 
given values of E and f. As our goal is to build a parameter estimation scheme 
based on linear or cubic spline approximations for the state, we first transform 
the domain n into a new domain where spline basis elements are more easily 
constructed. In this case it is natural to transform from Cartesian to polar 
coordinates so that n is the rectangle (e,R) ^ (0,2n); spline elements are easily 
defined on fi, whereas on one must handle the problem of defining such elements 

on the circular boundary of the annulus. 

We transform from Cartesian to polar coordinates in the usual way; appealing 
to Theorem 3.35, p. 63 of [26], we note that properties of this transformation 
guarantee the existence of a bounded operator T : H^(q) -»■ H (q) that is merely 
the composition map 

(Ty)(r,o) - y(4)(x^,X2)) 


^81) 


37 


where (j) : n ->■ i s the Cartesian to polar coordinate transformation. From [26], 
we also have that T ^ exists and is bounded. In polar coordinates, the state 
equation (73) becomes an equation on fi, 

a^(y,v) = <F,v>j^ (82) 

where a^(y,v) = // {TE)v(Ty) • v(Tv)do 

n 

and <F ,v> h // (TF) • (Tv)dfi, (dfi = rdrde). Henceforth we shall not 

distinguish between <,> and <,>,^ or a vector w and its transformed form Tw; 

however it should be kept in mind that in (73) both y and v are in hJ(«) while 

in (82) the transformed state Ty and associated vector Tv are in 

^0,per^^^ = THq(q). It is worthwhile to note that for smooth y€ Hg(fi) (so that 

we actually have y = 0 on 3fi), the. transformed y e hI (ii) satisfies 

u,per 

y(e,o) = y(R,G) = 0 and the periodic condition y(r,0) = y(r,2n), (r,e)€ fi . 

In addition, we remark that the gradient operator 7 in (82) is the polar 
coordinate transformation of the usual gradient in Cartesian coordinates; that 
is, in (82), 

w^cose - -p WgSino \ 

I (83) 

w sine + w.coso / 
r r 0 

for w€ H^(q). It is interesting to note that although v in (83) is the correct 
form of the gradient in polar coordinates, one often sees in the literature the 
operator 7 given by 

( 84 ) 





38 


In fact V and 7 may be used interchangeably in some circumstances (and indeed 
the latter greatly simplifies calculations), largely due to the fact that 


VW • VZ = 7W • 72 (85) 

and 

7 • 7W = V • 7W = AW (86) 

for w, z e H^(n). One must be careful however to use 7, and not 7, in cases 
involving an integration by parts (Green’s formula). 

We turn now to spline-based approximation for the state variable y in (82). 
For positive integers M and N we subdivide [e,R] and [0,2n] by defining partition 
points = e + (R-e)(i-l )/M, i = 1, .... M+1, and o'J = 2n(j-l)/N, 
j = 1, ..., N+1. We then construct a k-th order spline approximation y^’^^ for 
y using Galerkin techniques. That is, we find = span 

i = 1, .... M, , j = 1, ..., N, } where the basis elements v^’^. are tensor 
products of the k-th order spline elements, 


M,N 


* 1 J 


(87) 


M N 

both aj^ ^ and 6|^ j are standard k-th order B- spline elements that have been 
modified to match the Hq pg^(^^) boundary conditions. That is. 


«]^,i(c) = «k,i(R) = 0 , all i, 
and 

3|Jj(0) = 3kj(2ll), all j. 


I 
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We remark here that the dimension of depends on k as well as on M and N; 

in the case of linear splines (k = 1 ), = M - 1 and = N. 

To approximate the solution y to (82) for fixed k and Ee Q^, f e given, 
we look for a solution in to the Galerkin equations, 




( 88 ) 


for i = 1, .... and j = 1, n^. In the case of linear splines (k = 1), 

substitution of (87) into (88) yields the algebraic equations found in (15) of 
Section IV. There is no difficulty in establishing existence and uniqueness of 
solutions to (88) given the properties of the inner product matrices that appear 
in the algebraic equations derived from (88). 

Our goal is to construct a sequence of parameters (E*^’^, f*'^’*'*) 6 x 
and associated solutions ?”•'*) to (88) and detennine that 

JM,N, ^ ?«•") , y(E,f) in son,e sense; here (E.f) is 

an "optimal" parameter vector in the sense that it minimizes a distributed least 
squares fit-to-data criterion of the form 


J(F.f) = // I y(E,f) . d$2. 


where y^^ is the interpolated form of the pointwise data (transformed to 
satisfy homogeneous boundary conditions, as defined in (54)). If we use the 
form of the cost functional given in (57) we must strengthen the convergence 
results given in this section which can be done, but at the expense of additional 
technical tedium. Given pointwise data, we can use the J just defined if we 
interpolate this data. The minimization of the cost functional is performed 
over some constraint set Q > F. In the case of linear splines, we choose 
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Q X F “ Q-| later results will show, we do obtain the desired 

parameter and state variable convergence as M,N ®. If, instead, cubic splines 
are to be used, convergence results are obtained at the expense of the require- 
ment that Q and F be smoother sets. In particular, we need Q - and f = ^3 
where, for p = 2, 3, 

Qp = {E e cP(n) |E € Q^, |Elp < 1) 

Fp S {f e I f e F^, |f|p_^ < m) . 

Before considering parameter and state convergence, we establish the con- 
tinuous dependence of y, yj^’^ on the parameters E and f. 

Lemma 2: Let k = 1 or 3. 

(a) The mapping (E,f) -> y|^’^(E,f) : x f.| ->■ HQ(fi) is continuous in the 

I 2 X L^ topology on x f-j . 

(b) The mapping (E,f) y(E,f) : x Hj(fi) is continuous in the 

L^ X L 2 topology on x f^ . At (E,f) £ Q 2 x ^2’ ^ continuous in 
the L 2 X L 2 topology on x f^. 

Proof : We first consider the continuity of (E,f) -> y(E,f). Let E, E €. Q-j , 

f, f €. F^ and let y = y(E,f), y = y(E,f). Then, 

// Ev(y-y) • vvdi> = //(E-E)vy • vvdQ + // (Evy - Evy)-* vvdn 

■ h 2 

= If (E-E)vy • vvdn + a^(y,v) - a^(y,v) 

ii 

As ^ 

= // (E-E)vy • vvdfi + <F,v> - <F,v> 
n 



41 


for all v€ H^in) where we have used (82) in the last equality and 

F = f + V . (EvU), F = f + V . (EvU). In particular, if we let v = y - y, 
we obtain 

// F |v(y-y)| dfi = fj (E-E)vy • v(y-y)dn + <f-f, y-y> 

+ <v • {(E-E)VU}. y-i> . 

After an integration by parts and several applications of the Cauchy- Schwarz 
inequality, we have 

m|v(y-i)|^ < |(E-E)vy| |v{y-y)| + |f-f| |y.;| + j (E-E)7U| | v(y-i) | 

iK^ |(E-E)vy| + Kplf-fl + K^|7U|^ |E-E| (89) 

where here |*| denotes the L 2 (n) norm, !♦! denotes the L (n) norm and K„, K, 

L. cw 00 0 I 

are given in Lemma 1. Thus 

|y-y|l 1 c(|(E-E)7y| + |f-fj + |E-E1) (90) 

for some constant c > 0. For (E,f) e x we have, since y< H^, 

|y-y|f 1 c(|E-El^ |v^| + |f-f| + c|E-E|^) 

< c(K^|i-E|^Mf-f| + 3|E-E|J . (91) 

so that (E,f) y(E,f) is continuous in the x topology. If in addition 
(E,f)€ X SO that y € H (from the proof of Corollary 1), we obtain from 
the Sobolev Imbedding Theorem (see p. 97, [26]) 
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KyLic|y|3 

where is In H because (E,f) € x (again see the proof of 

Corollary 1). The continuity of (E,f) > y(E,f) at (E.f) in the L 2 x L 2 
topology is thus assured since (90) may now be rewritten as 

|y-yli 1 c(|E-E| |vi|^+ |f-f| + |E-i|). ( 92 ) 

For k = 1 or 3, the continuity of (E.f) ^ y5J’^(E,f) in the L 2 x L 2 
topology is apparent from (88). (See, for example, the algebraic eguations 
(15)-(20) in Section IV that replace (88). in the case of k = 1; we note that 

C|f 1 for i,p = 1, M-1 and j,f = 1, N, where 

C is a constant depending on M and N. A comparable statement may be made about 
the matrices involved in the algebraic equations derived for the case of k = 3.) 

Our later convergence arguments rely heavily on certain spline estimates 
that are simple modifications of those found on pp. 83-84 of [20a]. (In particular, 

our calculations must take into account the fact that the spline basis elements 

M 

**k,1 3nd on [0,2n] satisfy zero and periodic boundary conditions 

respectively. This causes no difficulty since we will only use those elements 
to approximate functions that satisfy such boundary condtions. For example, 
in the case of linear splines, {a^.} contains all the standard B-spline "hat 
function" basis elements {a'!'} except the "half-hats", Sl’J, If geHj(e,R), 

then the linear interpolating spline for g at knots (r^) is given by 
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PASS JS 

OF. POOR QUALfTY 


N •'vN 'vJ'i ' 

one element given by We thus net, for a ?:i-periodi'.; function g, 




“ y + (d^ + ] ) 3|,j • ) The spl ine estimates needed for 

later convergence proofs are given below. 

Lemma 3: Let y» e O h 1 (i») for k = 1 or 3. If denotes the 

*■ U j P vi I K 

canonical projection of L,(q) onto (in the L^iw) topology), then 

dM.N., 


'I'l 1 CQ(MN)"^''^^^|D‘'^^:;.i 

(93) 

^’^)| < C^(MN)‘''lD‘^'',j.! 

( 94 ) 


(95) 


Finally, we consider the existence of "optimal" parameters (L,T), 
detennination of a seqiionce { (E^’^,f^'’*^) 1 of approxinicit inn naramerers urid the 
convergence, in an appropriate sense, of to (i:,0 ao;! tiu? corresponding 

state variable convergence, of to ya',?). 

Le mma 4 : Let Q.j < F.j be compact in the LoCi-) > ^ P-'Oduct topology and 

let k = 1 or 3. Then there exists a solution A U, • F, to the 

! i 

M N 

problem of minimizing over ■ F, where 
J^’^E,f) = // iy'J;’^E,f) - yj^da 

5»": 

and y^|^’*^{E,f) is the solution to (88) associated witli parametet'S (E,f). 

Proof : The result follows immediately from Lemnia 2, which guarantees the 

M N 

continuity of over the compact set Q.j >= F^ . 
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The sense in which approximates (E.f), and in turn 

approximates y(E,f), is given in our final two results. 

T h eorem 2 : Let { } be a sequence of parameters in Q, x f and let 

M N ' * 

y|^* be the solution to (88) corresponding to E*^’^,f*^’*^. 

(a) If there exists (E,f)€ x such that (E^’^f^’^j ^ 

X L2 topology, then y'J’^ y(E,f) in Hj(fi) as M.N ■> «. 

(b) Suppose further that (E^’^ f^’^j e x f^ , (E,f)e x and 

that ^ (E,f) in the L2 x L2 topology. Then ylj'*^'' y(E,f) 

in hJ('^) as M.N -> ... if (E,f) e Q3 v F3, then y^’^ 

HgCn) as M.N ■>■ «. 


Proof: For parts (a) and (b) it suffices to show that 

- y(E,f) 1^ -> 0 as M.N -*■ «>, 


uniformly in (E.f) (in the desired parameter set), because we may then use this 
result to demonstrate that 




* - y(E.f)|, 

where both terms converges to zero as M.N ^ ». Here we have used the appropriate 
continuity result from Lemma 2 to argue the convergence of the second term above. 
Let (E,f)e X F^. For k = 1 or 3 we subtract (88) from (82) (with 

“MJ « "J.per*")) obUiu 
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aE(y(E,f) - = 0 

for i = 1, ..., M, , j = 1, N,. That is, is the projection of y€ 

M N . 1 

onto Sj^’ in the a^(*,*) topology (equivalent to the Hq topology). We thus 

obtain, letting y £ y(E,f), y”’'^ £ y5J’^(E,f), 

>n|v(y - yJJ’’^)!^ < a^(y - y - y”’^) 

< a^(y - y - p”‘%) 

< iniv(y - p[J’'^y)i^ 

; rac{|D^(y - f |D.(y - P^’%)1^1 (96) 

M N , . 

where. P|^’ was defined in Lemma 3 and (83) was used to establish the last 
inequality. 

We consider part (a). For k = 1 and (E,f)€ q.^ fact 

that y€ H^(q) (Corollary 1) and estimates (94), (95) and (75) to claim that the 
bound in (96) becomes 

l^(y - y'J'’'^)|^ < 2C^^(MN)’^ ^ |D^y|^ 

1 C(MN)'^(|f| + |E|i + |y|)2 

for some constant C. From the definitions of Q^, and the a priori bound (77) 
on y (uniform in (E,f)) we thus have that 

|v(y - y1’’^)|^ - C(HN)"^ 

for some constant C, where C is determined independent of M, N, E and f. 
Therefore, [vT ^(y - y*^’*^)|^ = 0((MN)'^) so that (writing y, y^*^ for T"V» 
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^ y in Hj(fi) as M.N -v », uniformly in (E.f) € x Hence 
part (a) obtains. 

To prove part (b), we repeat the above arguments for the case of k = 1 
and similarly obtain the convergence of y^*^ to y in HQ(fi), for fixed 
(E.fOC" Q 2 X To consider the case of cubic spline approximations (k = 3), 
we note that (E,f)£ x F^ implies that ye H^(n) (Theorem 2) and that, 
following (96) and the arguments for part (a), 

lv(y - y3^)|^l2C^(MN)‘^^ |0^i^ 

1 C(MN)’^(|f |2 + IEI 3 + ly|)^ 

using (76) to obtain the second bound. Using (77) and the definitions of 03 * ^ 3 * 
we are able to claim that 

|v(y - < c(HN)'^ 

so that y^* ^ > y 1o hI(o} unifonnly in (E.f) € Q.^ . The proof of part (b) 
is then complete. 


Theorem 3 : (a) Let x f^ be compact in the x L 2 topology. For each 

M,N there exists a solution (E*^*^,f^’*'*) to the problem of minimizing over 
X F-|. In addition, there exists (E,f)£ x f-j and a subsequence (relabelling) 

{(i^j’^j, of {(E^*'^,f^’'^)} such that (§^J ,f^J (E,f) in x L 2 , 

^ y(E,f) in HQ(fi), and (E,f) is a solution to the original 
problem of minimizing J over F^. 
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(b) Let k = 1 or 3. Let xf be compact in the L. • L, topology 

and (E ,r ) a solution to the problem of minimizing over ^ 

There exists (E,f)C Q.jxy^ and a subsequence {(E such that 

-Mj , N j -Mj » Nj - - 

‘■f t ^ (E.f) in L 2 X L 2 * If in addition (E,f)€ x then 
Mj.Nj - Mj.Ni -M-.N. - . , 

^ -*■ y(E.f) in Hq(s^i) and (E,f) is a solution to the original 

parameter estimation problem of minimizing J over x if (E,f) € x f^ 
the same statement may be made in the case of k = 3 (cubic spline approximations). 

Proof; It is easy to argue, as in the proof of Lemma 2(a), that 
(E.f) y^’^(E,f) is continuous in the L^^ x L 2 topology, so that is 

continuous over the compact set x f, . Thus, for each M.N, there exists a 
solution (E ’ ,f ’ ) to the approximate parameter estimation problem. From 
the compactness of x f^ , there exists a subsequence and 


• I ■ ' — 1 ’ I 

an element (E.f ) e 0^ such that > {E.f) in L L 

from Theorem 2, . y(E,f) in hI(JS). Finally” 


and. 


w- Mi.N. M.,N. 

J(E,f) = lim J J J(E ^ J f J J\ 

J 0 


lim 


V"o- ’ 


Mi,Ni 

-^(E,f) 


= J(E.f) 


for any (E,f)€ x f^ , where we have used Theorem 2(a) and the fact that 
(E J J,f J J) i; 


M ■ ,N • 

is a minimizer for ^ 


over x f^. Therefore, (E,f) 
is a solution to the original parameter estimation problem. 
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To prove part (b), we duplicate the above arguments for both linear and 
cubic spline approximations, applying Theorem 2(b) and the L 2 topology on 
Qg (Q3) throughout. 


Before our theoretical justification of the spline-based estimation scheme 
is complete, we must indicate one final result that focuses on the means by which 
we solve the (M,N)-th approximate parameter estimation problem.. It is implicitly 
assumed, in the previous calculations that we are able to compute , a 

solution to the problem of minimizing over a suitable parameter set Q x F. 

In fact, the sets Q x y considered here are infinite-dimensional so that we 
must perform the minimization over a function space. To implement the 
approximate estimation scheme, we actually p>arametrize E and f and minimize 

over a finite dimensional set x f^"^. To see how (e*“*»*'*»^ »J) 

( a solution to this problem) approximates (E,f), we are forced to consider the 
limit of such a sequence as M, N, I, J -»• «. That we actually obtain the theoretical 
convergence of such a sequence (or subsequence) is an easy extension of the ideas 
in [15a]; these techniques have already been successfully implemented and are 
described, for the case of the antenna problem, fn section IV. 
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FIGURE 2. Side View of Maypole (Hoop/Column) Antenna 
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Figure 9: Example 3.5 
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